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Abstract 

Exclusive electroproduction of mesons has been measured using the ZEUS 
detector at HERA in two ranges, 0.25 < < o.85 GeV^ and 3 < < 30 GeV^. 
The low-Q^ data span the range 20 < < 90 GeV; the high-Q^ data cover the 
40 < < 120 GeV interval. Both samples extend up to four-momentum transfers 
of \t\ = 0.6 GeV^. The distribution in the azimuthal angle between the positron 
scattering plane and the production plane shows a small but significant violation 
of s-channel helicity conservation, corresponding to the production of longitudinally 
polarised (i.e. helicity zero) p^ mesons from transverse photons. Measurements of 
the 15 combinations of spin-density matrix elements which completely define the 
angular distributions are presented and discussed. 
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1 Introduction 



Exclusive production of vector mesons by real and virtual photons, 7p — > Vp and 7*]? — >■ 
Vp {V = p°, uj, (j), J/ip, . . .), has been studied extensively, both in fixed-target experiments 
and at HERA |^. For photon-proton (7*p) centre-of-mass energies > 10 GeV, the 
production of light vector mesons (p°, u, (p) at low photon virtuality, < 1 GeV^, exhibits 
features typical of soft diffractive processes: a weak dependence on centre-of-mass energy 
and a differential cross section that, at low \t\ values, falls exponentially with —t, where 
t is the square of the four-momentum exchanged between the photon and the proton. 
These features are consistent with the expectations of the Vector Meson Dominance model 
(VMD) 0] according to which the photon is considered to fiuctuate into a vector meson 
which then scatters elastically from the proton. 

The low-energy data indicate that the amplitude for the photon/ vector- meson transition is 
predominantly s-channel helicity conserving, i.e. the helicity of the vector meson is equal 
to that of the photon when the spin-quantisation axis is chosen along the direction of the 
meson momentum in the 'j*p centre-of-mass system. The exchange of an object in the 
t-channel with P = (— 1)"^, e.g. = 0~^, 1~, 2"*", 3~, . . .(natural parity exchange in the t 
channel, NPE) , dominates such diffractive processes Q . However, small helicity-single- 
flip and helicity-double-flip contributions to the amplitude have been reported in tt+tt" 
photoproduction in the p° mass region at W < 4 GeV [H. Helicity-single-flip amplitudes 
have also been observed in p° electroproduction for 1.3 < W < 2.8 GeV and 0.3 < 
Q'^ < 1.4 GeV^ p. A helicity-single-flip contribution of (14 ± 8)% was measured in p° 
muopro duct ion at W = 17 GeV 

Data from HERA have shown that when the reaction involves a large scale, such as the 
charm quark mass in J/ ip photoproduction [§] or high in p° electroproduction |10 



the cross section rises more steeply with W and the t distribution is flatter than in the case 
where no such scale is present, such as p° photoproduction |Tl], [T2|, |T3[. Such results can 
be understood in terms of models based on perturbative QCD (pQCD). In these models, 
the photon fluctuates into a qq pair and the interaction of the pair with the proton is 



mediated by the exchange of two gluons in a colour singlet state fT^ . The resulting cross 



section is proportional to the square of the gluon density in the proton. 

The cross section for the exclusive production of vector mesons from virtual photons has 
contributions from both transverse and longitudinal photons. These contributions are 
expected to have different W and dependences |l|, |1^, and their extraction would 
provide an important test of the understanding of the production process. The angular 
distributions of the decay products of the vector meson, in principle, yield the separate 



contributions of transverse and longitudinal photons. Previous studies at HERA 0, [IT 



|T^, [T^, assumed conservation of helicity in the s-channel amplitudes (s-channel 

helicity conservation, SCHC), in which case the ratio of longitudinal to transverse photon 
cross sections, i?, is given in terms of the ratio of the number of helicity (i.e. lon- 
gitudinally polarised) to helicity ±1 (i.e. transversely polarised) p° mesons produced. 
Sufficient data are now available to test the validity of SCHC at HERA by measuring the 
full set of matrix elements which completely determine the decay angular distributions. 
This paper reports the determination of such a set of matrix elements for exclusive p° 
electroproduction at HERA, ep — > ep^p (p° tt+tt"). Two kinematic ranges are studied. 
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0.25 < g2 < 0.85 GeV2, 20 < < 90 GeV (referred to as the "BPC" sample) and 
3 < < 30 GeV^, AO < W < 120 GeV (the "DIS" sample). The two samples are those 
already used in the previous ZEUS paper |P, where SCHC was assumed. A similar analysis 
has recently been presented by the HI Collaboration in the range 1 < < 60 GeV^ |T0|. 



The data are used to determine R without relying on the SCHC hypothesis, to evaluate 
the size of the helicity-single-fiip and helicity-double-fiip amplitudes and to test the extent 
to which natural parity exchange dominates in the t channel. The results are compared 
to recent calculations [|19| , pO| , of the amplitudes for exclusive p° meson production in 
the framework of pQCD models of the type described above, i.e. assuming the t-channel 
exchange of a gluon pair. 

Under the assumption of NPE, there are five independent combinations of photon and me- 
son helicity states, yielding two helicity-conserving amplitudes, two single-flip amplitudes 
and one double-flip amplitude for p° electroproduction (see, for example, []19[)- The cal- 
culations indicate that, for photon virtualities exceeding the hadronic mass scale of about 
1 GeV^, the single-flip amplitude for producing longitudinally polarised vector mesons 
from transverse photons is significant. This amplitude is of higher twist and vanishes 
in the non-relativistic limit of the p° meson wave function; in addition, it is propor- 
tional to \f\t\, while the double-flip amplitude grows linearly with \t\. These amplitudes 
are rendered experimentally accessible via their interference with the dominant helicity- 
conserving amplitudes. Conversely, the amplitude for producing transversely polarised p° 
mesons from longitudinal photons is expected to be negligible There are also predic- 
tions of non-vanishing single and double-flip contributions in the non-perturbative region 
of small photon virtualities, < 1 GeV^ ^ . VMD-based models |2^ also predict 



SCHC breaking on the basis of the analogy between photon-hadron and hadron-hadron 
interactions. 



2 Experimental set-up 

The data were collected at the ep collider HERA in 1995 using the ZEUS detector. In this 
period HERA operated at a proton energy of 820 GeV and a positron energy of 27.5 GeV, 
giving a total centre-of-mass energy a/s — 300 GeV. A detailed description of the ZEUS 
detector can be found elsewhere 0, The main components used in this analysis are 
briefly described below. 

The high-resolution uranium-scintillator calorimeter, CAL, consists of three parts: for- 
ward [| (FCAL), barrel (BCAL) and rear (RCAL) calorimeters. Each part is subdivided 
transversely into towers that are segmented longitudinally into one electromagnetic sec- 
tion and one (RCAL) or two (FCAL, BCAL) hadronic sections. 

^Throughout this paper the standard ZEUS right-handed coordinate system is used: the Z-axis points 
in the direction of the proton beam momentum (referred to as the forward direction) and the horizontal 
X-axis points towards the centre of HERA. The nominal interaction point is a,t X = Y = Z = Q. 
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Charged-particle tracks are reconstructed and their momenta measured using the central 
tracking detector (CTD). The CTD is a cyhndrical drift chamber operated in a magnetic 
field of 1.43 T produced by a superconducting solenoid. It consists of 72 cylindrical layers, 
organised in 9 superlayers covering the polar angular region 15° < 9 < 164°. 

The trajectory of positrons scattered at small angles with respect to the beam direction 
is determined from the beam pipe calorimeter (BPC) and the small-angle rear tracking 
detector (SRTD). The BPC is an electromagnetic sampling calorimeter located aX Z = 
—294 cm. The SRTD is attached to the front face of the RCAL. It consists of two planes 
of scintillator strips, 1 cm wide and 0.5 cm thick, arranged in orthogonal orientations and 
read out via optical fibres and photomultiplier tubes. It covers the region of 68 x 68 cm^ 
in X and Y with the exclusion of a 10 x 20 cm^ hole at the centre for the beam pipe. 



3 Kinematics and decay angular distributions 

The following kinematic variables, some of which are indicated in Fig. ^ are used to 
describe exclusive p° production: 

• The four-momenta of the incident positron (A;), scattered positron (/c'), incident 
proton (P), scattered proton (P') and virtual photon (g); 

• Q"^ = —q^ = —(k — k'Y, the negative four-momentum squared of the virtual photon; 

• W"^ = {q + Py, the squared centre-of-mass energy of the photon-proton system; 

• y = {P ■ q)/{P ■ k), the fraction of the positron energy transferred to the proton in 
its rest frame; 

• Mtttt, the invariant mass of the two decay pions; 

• t = {P — P')^, the squared four-momentum transfer at the proton vertex. 

The kinematic variables were reconstructed using the so-called "constrained" method 
25|1 , which uses the momenta of the decay particles measured in the CTD and the polar 



and azimuthal angles of the scattered positron determined with the BPC (for the BPC 
sample) or with the CAL and the SRTD (for the DIS sample). 

The exclusive electroproduction and decay of p° mesons is described, at given values of 
W, Q^, M^TT and t, by three angles: - the angle between the p° production plane and 
the positron scattering plane in the 7*p centre-of-mass frame (see Fig. 0); 6h and (ph ^ the 
polar and azimuthal angles of the positively charged decay pion in the s-channel helicity 
frame, in which the spin-quantisation axis is defined as the direction opposite to the mo- 
mentum of the final-state proton in the p° rest frame. In both the rest frame and the 
7*p centre-of-mass system, (p^ is the angle between the decay plane and the p° production 
plane. The angular distribution as a function of these three angles, W {cos 9h,(ph,^h), is 
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parameterised by the p° spin-density matrix elements, p"^, where i,k = —1,0, 1 and by 



convention a=0, 1,2,4,5,6 for an unpolarised charged-lepton beam The superscript 
denotes the decomposition of the spin-density matrix into contributions from the following 
photon polarisation states: unpolarised transverse photons (0); linearly polarised trans- 
verse photons (1,2); longitudinally polarised photons (4); and from the interference of the 
longitudinal and transverse amplitudes (5,6). 

The decay angular distribution can be expressed in terms of combinations, and r^^, of 
the density matrix elements: 



ik 



Pik + (^Rpik 

1 + eR ' 

Pik 

1 + eR' 

"RpTk 

1 + eR' 



a 



a 



1,2 



5,6, 



(1) 
(2) 



where e is the ratio of the longitudinal to transverse photon fluxes and R = ^/aj^^, 
with (t2 ^ and ctj. ^ the cross sections for exclusive p° production from longitudinal and 
transverse virtual photons, respectively. In the kinematic range of this analysis, the value 
of e varies only between 0.98 and 1.0; hence p^^, and pf^ cannot be distinguished. 

The hermitian nature of the spin-density matrix and the requirement of parity conserva- 
tion reduces the number of independent combinations to 15, in terms of which the angular 
distribution can be written as 



3 

W {cos Oh, (l)h,^h) = j- 

4:71 



2^ 



"oo) + ^(3rSo'-l)cos^^. 



-\/2 Rejr^Q} sin 26 h cos (f)h — r^l^ sin^ 6h cos 20/j 

-e cos 2^h{rl^ sin^ 6h + "^^s^ 6^ — V2 Re{r|Q} sin 26 ^ cos < 

-rl^^sin^ 6h cos 20,,) 
-e sin 2^h{V2 Im{r^Q} sin 26h sin 0^ + Im{r^_]^} sin^ 6h sin 2 



+ J2e{l + e) cos ^hirti sin^ 6h + r^o cos^ 6h 



-v/2 Re{r^o}sin2^;, cos (j)h — rl_^ sin^ 6h cos 



+ y^2e(l + e) sin $,,( lra{r%} sin 26h sin (ph 
+lm{rl_i} sin^ 6h sin 20, 



(3) 



The 15 coefficients r°^, are related to various combinations of the helicity amplitudes 
T\pX^, where Xp and are the hehcities of the p° meson and of the photon, respectively. 
Table |] shows the relations between the coefficients r°^, rfj^ and the helicity amplitudes. 

Under the SCHC assumption, the angular distribution for the decay of the p° meson de- 
pends on only two angles, 6^ and iph = (ph — ^h, and is characterised by three independent 
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parameters, rQQ,rl_^ and Re{rfQ}: 



W{coseh,iph) 



-I- 



1 



(3rOo'-l) cos' 9, 



+ erl_^ sin' 6 h cos 211; h - 2y^e(l + e)Re{r^o} sin 26'/^ cos?/;^]. (4) 

The SCHC hypothesis also imphes that r\_-^ = —lm{rf_^} and Re{r^Q} = — Im{r^Q}. In 
this case, the ratio R can be determined from the polar angular distribution alone: 



R 



„04 
00 



P 1 _ 7-04 ■ 



(5) 



The combined assumptions of SCHC and NPE further reduce the number of independent 
parameters to two, since the polar and azimuthal angular distributions are then related 
via: 



1 - 4' - 2rl_i = 0. 



If the SCHC requirement is relaxed, then 



R=- 



„04 
' 00 



e 1 - (rgo^ - A2) 



(6) 
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with 



where Nl 
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|TiiP + IT, 
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+ |Tc 



oil 



(8) 

The quantity A can be 



approximated as: 
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oil 



ITf 



oil 



00 



y/Nr + eNL 



TiiP + lTc 



00 



(9) 



where the amplitude |Tio| (for the production of transversely polarised p° mesons from 
longitudinal photons) is neglected with respect to |Toi| and the amplitudes Tqo and Tqi 
are taken to be in phase; in addition e = 1 is assumed. Single and double flip amplitudes 
have been assumed to be small with respect to Tqo and Tn. 



Finally, the assumption of NPE alone leads to the sum rule: 

+ 2r°li - 2ri\ - 2rJ_i = 0. 



1 - 

'00 



(10) 



4 Event selection 



The event selection is similar to that of the previous analysis 0; it is summarised briefly 
in the following. 

The online event selection required an energy deposit in the BPC or an electron candidate 
in the CAL, along with the detection of two tracks in the CTD. 

In the offline selection the following further requirements were imposed: 
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The energy of the scattered positron was required to be greater than 20 GeV if 
measured in the BPC and greater than 5 GeV if measured in the CAL. 



The requirement E — pz > 45 GeV was imposed, where E — pz = Y^ii^i — pzj 
and the summation is over the energies and longitudinal momenta of the final-state 
positron and pions. This cut, apphed in the DIS analysis, excluded events with high 
energy photons radiated in the initial state. 

The Z coordinate of the interaction vertex was required to be within ±50 cm of the 
nominal interaction point. 

In addition to the scattered positron, the presence of exactly two oppositely charged 
tracks was required, each associated with the reconstructed vertex, and each with 
pseudorapidity |?7| less than 1.75 and transverse momentum greater than 150 MeV; 
pseudorapidity is defined as rj = — In (tan (^^/2)). These cuts excluded regions of 
low reconstruction efficiency and poor momentum resolution in the CTD. 

Events with any energy deposit larger than 300 MeV in the CAL and not associated 
with either the pion tracks or the positron were rejected 



In addition, the following requirements were applied to select kinematic regions of high 
acceptance: 



• The BPC analysis was limited to the region 0.25 < Q'^ < 0.85 GeV^ and 20 < ly < 
90 GeV. 

• The DIS analysis was restricted to the kinematic region 3 < < 30 GeV^ and 
40 < < 120 GeV. 

• Only events in the tt+tt" mass interval 0.6 < M.^„ < 1.0 GeV and with \t\ < 0.6 GeV^ 
were taken. The mass interval is slightly narrower than that used previously 0, in 
order to reduce the effect of the background from non-resonant tt+tt^ production. 
In the selected M^^.,^ range the resonant contribution is ~ 98% for the BPC sample 
and ~ 100% for the DIS sample. 



The above selections yielded 5271 events in the BPC sample and 2510 events in the DIS 
sample. 



5 Monte Carlo simulation 



The relevant Monte Carlo generators have been described in detail previously Here 
their main aspects are summarised. 



In the BPC analysis, a dedicated Monte Carlo generator based on the JETSET7.4 P9 
package was used. The simulation of exclusive p° production was based on the VMD 
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model and Regge phenomenology. The effective Q^, W and t dependences of the cross 
section were parameterised to reproduce the measurements 0. 

In the DIS analysis, a program interfaced to HERACLES4.4 |3l|] was used. Again, the 
effective Q^, W and t dependences of the cross section were parameterised to reproduce 
the results 0. 

In both cases the decay angular distributions were generated uniformly and the Monte 
Carlo events were then iteratively reweighted with Eq. (|^) using the results of the present 
analysis for the 15 combinations of matrix elements r°^, r^. 

The generated events were processed through the same chain of selection and reconstruc- 
tion procedures as the data, thus accounting for trigger as well as detector efficiencies and 
smearing effects. For both analyses, the number of simulated events after reconstruction 
was about a factor of seven greater than the number of reconstructed data events. Fig- 
ure 1^ shows the acceptances (defined as in Ref. as a function of cos 9h, (ph and 
both for the BPC and the DIS data. The different normalisations of the acceptances for 
the BPC and DIS samples mainly reflect the different coverage for the azimuthal angle of 
the scattered electron, the BPC covering only a small fraction of 2tt. 



6 Results 

Figures ^ and ^ show the reconstructed decay angular distributions for the data (full 
symbols) and Monte Carlo events (solid histograms) for the BPC and for the DIS samples, 
respectively. The Monte Carlo simulation reproduces the measured distributions well. The 
dashed histograms show the same Monte Carlo events reweighted to respect SCHC. The 
violation of SCHC can be most clearly seen in the distributions, particularly for the 
DIS sample of Fig. ^ In the following we focus on how this violation manifests itself in 
terms of the 15 combinations of spin-density matrix elements, r°^, rj^. 

The combinations r°^, were obtained by minimising the difference between the three- 
dimensional {cos 9h, (f>h,^h) angular distributions of the data and those of the simulated 
events. A binned likelihood method was used, assuming a Poisson distribution for the 
number of events in each bin. The number of bins in cosOh, 4>h and $/j was 8x8x8. 

Table | summarises the results for the 15 coefficients t^^, rff,. The BPC results cover 
the kinematic region 0.25 < < 0.85 GeV^ {{Q^) = 0.41 GeV^), 20 < 1¥ < 90 GeV 
({W) = 45 GeV), 0.6 < M^^ < 1.0 GeV and \t\ < 0.6 GeV^ ((|t|) = 0.14 GeV^). The 
DIS results cover the region 3 < < 30 GeV^ [{Q^) = 6.3 GeV^), 40 < < 120 GeV 
({W) = 73 GeV), 0.6 < M^^ < 1.0 GeV and |t| < 0.6 GeV^ ((|t|) = 0.17 GeV^). 
Tables ^, || and ^ give the BPC results in bins of W, and t. The limited statistics and 
poorer t resolution do not allow such fine binning for the DIS sample. Tables ^ and |^give 
the full correlation matrix for the BPC and the DIS results, respectively. 

Figure ^ shows the combinations of matrix elements for the BPC and the DIS data as 
a function of Q^. The coefficient r^g represents the probability that the p° meson be 
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produced in the helicity zero state, i.e. with longitudinal polarisation; it grows with Q^, a 
result qualitatively in agreement with SCHC, since the contribution of longitudinal pho- 
tons increases with increasing . Less pronounced variations with are also observed 
for r\_Y and Im{r^_^}, which are related to Tqq under the combined SCHC and NPE 
assumptions (see Sect. |^). 

Figure ^ shows the coefficients r\^, as functions of t for the BPC sample alone. The 
data suggest a growth of with \t\. 

Figures |, and | compare the results to the expectations of SCHC (dashed lines) for the 
ten elements for which SCHC makes explicit predictions. Again deviations from SCHC 
are observed for both the BPC data (RejrJo}, '^i-d ''^ii) I^^l^io}) ''^oo) ^^^1 the DIS data 
(rgg). A comparison of the results with the SCHC hypothesis yields x^/ndf = 81.1/12 
for the BPC sample and 'x^ /ndf = 79.2/12 for the DIS sample. Figures ^ and § also 
show the recent HI results obtained in the kinematic region 1 < < 60 GeV^, 
30 < W < 140 GeV and \t\ < 0.5 GeV^; they are in excellent agreement with the present 
data. The results of the model calculations |T^, ^ are also shown in Figs. ^ and ^. 



The present results are not corrected for the contribution of non-resonant tt+tt^ produc- 
tion and its interference with resonant tt+tt^ production via the so-called Soding 
mechanism ||32|. Hence, strictly speaking, they only apply to the reaction ep evr+vr"^ 
with M^jr in the quoted range. While the relative magnitude of the non-resonant contri- 
bution is small, the interference term, which changes sign at the nominal p mass value, 
may significantly affect the angular distributions. In order to assess the sensitivity of 
the data to changes in the selected M^^r range - and hence changes in the relative con- 
tributions of the non-resonant and the interference terms - the events with M^-tt < Mp 
and Mtttt > Mp were analysed separately. For the DIS sample, where the Soding effect 



is small [y, |I0[, the results thus obtained are consistent with those found for the whole 
Mtttt range, suggesting that the observed breaking of SCHC is indeed a feature of the 
reaction ep ep^p. The BPC data, however, exhibit some dependence on Mtttt, as shown 
in Table |. In particular, Tqq decreases for the higher M^jr selection, a result qualitatively 
consistent with expectations based on the interference between resonant and non-resonant 



TT+TT" production [Q. The observed variations of rl_^ and Im{r^_j^} reflect that of Tqq 
under the combined SCHC and NPE assumptions (see Sect. ^. Appreciable M.^^, depen- 
dences are also observed for Tqq, Re{rfQ} and Imjr^Q}, which suggests an influence of the 
non-resonant tt+tt^ production on these combinations of matrix elements. 

The results are also not corrected for the contribution from the proton-dissociative reac- 
tion, ep — > ep^N , where is a state of mass (Mat < 4 GeV) which escapes undetected 
in the proton beam direction. The size of this background was estimated to be (23±8)% 
and (24^5)% for the BPC and the DIS samples, respectively. The W, Q^, cosOh, (ph and 
$ft distributions for the proton-dissociative events (tagged by activity in the FCAL 0) 
were found to be consistent with those for the exclusive events, as expected on the basis 
of the factorisation of the diffractive vertices [Q. It has therefore been assumed that the 



proton-dissociative contribution does not affect the decay angular distributions. 
Finally, any influence of radiative corrections on the angular distributions has been ne- 
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glected. 



7 Systematic errors 

The systematic uncertainties were obtained by modifying the analysis procedures as hsted 
below. 



(a) Sensitivity to the event selection criteria: 

• The Mtttt range was restricted to 0.7 < M-^^^ < 0.84 GeV for the BPC sample 
and 0.65 < M^^ < 0.9 GeV for the DIS sample. 

• The minimum positron energy as measured in the BPC was increased to 
23 GeV. 

• The minimum value of E — Pz was raised to 48 GeV. 

• The minimum track transverse momentum was increased from 150 MeV to 
200 MeV for the BPC analysis and to 300 MeV for the DIS analysis. 

• The pseudorapidity range was restricted to |?7| < 1.7 for the BPC data and 
|?7| < 1.5 for the DIS data. 

• The t range was restricted to 0.05 < \t\ < 0.5 GeV^. 

(b) Sensitivity to the Monte Carlo simulation: 

• The W distribution was reweighted by a factor W^, with k varying between 
-0.1 and 0.1. 

• The distribution was reweig hted by a factor 1/{1 + Q^/M^)'', where k varies 
between —0.2 and 0.2 and Mp is the p° meson mass. 

• The t slope was varied by ±1 GeV^^ for the BPC data and by ±0.5 GeV^^ for 
the DIS data. 



Sensitivity to the fitting procedure: the method of moments was used instead of 
the maximum likelihood method. Moments of the observed three-dimensional dis- 
tribution in the angles {cos 9h,4>h and ^h) were calculated, i.e. the distribution was 
projected on appropriate functions of cos 9h, 4>h and The same moments were 
evaluated for the reconstructed angular distribution from the Monte Carlo simula- 
tion and the difference between the moments in the data and those in the Monte 
Carlo was minimised by adjusting the values of the coefficients r°^, r°j^ used in the 
simulation. This method is a variation of that presented in Appendix C of Ref. PU 



(d) Sensitivity to the binning: the standard number of bins in cos 9h, <t>h and 8 x i 
was changed to 6 x 8 x 8 bins and, in the BPC case, also to 10 x 10 x 10. 
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The dominant effects come from the sensitivity to the binning, the modification of the 
Mtttt range and the fitting method (maximum hkehhood vs method of moments). As an 
example, for Tqq, the systematic uncertainty due to the binning is 26% (14%), and that 
due to the M^^ range is 15% (11%) for the BPC (DIS) analysis. 



8 Discussion 



In this section the r?^, rf^^ combinations which exhibit a deviation from SCHC are individ- 
ually examined. The implications of the measured violation of SCHC for R = aj^^/aj'^ 
are also discussed, as is the extent to which NPE dominates the present data. 

Deviations from the null values expected in the case of SCHC are observed for Relr^'g}, 
r?li, r}i, Im{r?o} and r^o- 



1. Re{rO^} = 0.034 ± 0.006 (stat.) ± 0.009 (syst.) for the BPC sample. The DIS 
sample gives Re{r°o} = 0.013 ± 0.010 (stat.) ± 0.022 (syst.), which is consistent 
with the BPC result but does not significantly deviate from zero. The HI DIS 
result, Re{r°^} = 0.011 ± 0.012 (stat.) (syst.), is consistent with zero [|10|. 



Non-zero values of Rejr^Q} have also been observed in low-energy photoproduc- 
tion 0, as well as in electroproduction and muoproduction experiments. 
Reference quotes values ranging from 0.01 to 0.16 for < < 1 GeV^ and 
1.7 <W < 2.8 GeV. At higher energy, 12.5 < W < 16 GeV, and for < 3 GeV^, 
the CHIO Collaboration found Re{r%} = 0.07 ± 0.03 0. 

This coefficient is sensitive (cf. Table |l]) to the interference of helicity-conserving 
and helicity-single-flip amplitudes, as well as to the interference of helicity-single- 
flip and helicity-double-flip amplitudes. Both single-flip amplitudes are involved 
here, i.e. that for the production of longitudinally polarised p° mesons from trans- 
verse photons and that for the production of transversely polarised mesons from 
longitudinal photons. 

The data are consistent with the predictions of |2y, ^ at low Q"^. Some disagreement 
with the predictions of [jl9|, pTl is visible at high Q^. 



2. r?ii = -0.040±0.009 (stat.) ±0.019 (syst.) for the BPC sample. No deviation from 
zero is observed for the DIS sample, r^l^ = 0.000 ± 0.011 (stat.) ± 0.008 (syst.). 
The HI DIS result, r^l^ = —0.010 ± 0.013 (stat.)^QQg3 (syst.), is consistent with 
zero [|T^]. 

Negative values for r^ti, consistent with the BPC result, though only marginally 
different from zero, were also measured in the range 1.7 < W < 2 GeV, < 
1 GeV^ §. A value of r°li = -0.07±0.03 was reported for 12.5 < W < 16 GeV [0. 
Negative values of p^-i were measured in photoproduction 

This term is proportional to the square of the helicity-single-flip amplitude for the 
production of transversely polarised p° mesons from longitudinal photons, Tio, and 
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is also sensitive to the interference of helicity-conserving and helicity-double-flip 
amplitudes. 

For this combination of matrix elements, as for Rejr^'o} discussed above, the relative 
contribution from longitudinal photons is an order of magnitude larger for the DIS 
sample than for the BPC sample, owing to the dependence of R (see Eq. p. 

The pQCD based models mil, |2l| reproduce the DIS data for rf_^. At low 
some disagreement between the model |21| and the data is observed. 



3. = -0.039 ±0.009 (stat.) ± 0.021 (syst.) for the BPC sample. Here again, no sig- 
nificant effect is visible in the DIS data, = -0.006 ±0.012 (stat.) ± 0.026 (syst.). 
No evidence for a deviation from zero of was found by the fixed-target experi- 
ments or by HI; the latter measured = —0.002 ± 0.034 (stat.) ± 0.006 (syst.) in 
DIS 0. 



This coefficient is sensitive to the interference between the helicity-conserving and 
the helicity-double-fiip amplitudes. 

Here also the models [|19|, |l| reproduce the DIS data. At low disagreements 
similar to those observed for r^ti between the model |^1[] and the data are visible. 

4. MrfJ = 0.039 ± 0.008 (stat.) ± 0.012 (syst.) (BPC). The DIS result, MrfJ = 
0.008 ±0.014 (stat.) ± 0.031 (syst.), is consistent with zero, although also consistent 
with the BPC result. The HI DIS result, Imfr^J = 0.023±0.016 (stat.)+|;;!|j9 (syst.), 
is also consistent with zero |]TU|. No significant deviation of Im{r^o} from zero was 



found in the earlier experiments. 

This coefficient is sensitive to the interference of helicity-conserving and helicity- 
single-flip amplitudes as well as to the interference of single-flip and double-flip am- 
plitudes. The single-flip amplitude involved is that corresponding to the production 
of longitudinally polarised mesons from transverse photons. 

The pQCD based models [1T9| , |l|] give a satisfactory description of the data. 



5. r^o = 0-051 ± 0.010 (stat.) ± 0.018 (syst.) (BPC) and r^o = 0-095 ± 0.019 (stat.) ± 
0.024 (syst.) (DIS). The HI experiment found r^^ = 0.093±0.024 (stat.)+J;[5;o (syst.) 
in DIS [|T^, in good agreement with the present ZEUS DIS result. None of the flxed- 
target experiments measured signiflcant deviations from zero for Tqq. 

This coefficient is proportional to the interference between the helicity-non-flip am- 
plitude for longitudinal photons Tqo and the helicity-single-flip amplitude Tqi for 
the production of longitudinally polarised p° mesons from transverse photons. The 
BPC data suggest that Tqq increases with increasing \t\. 

The deviation of r^Q from zero is directly related to the failure of SCHC to describe 
the $h distribution (cf. Figs. ^ and |^), as is apparent by integrating Eq. (^) over 
cos9h and (ph- 



Wi^f,) oc [l-ecos2$;,(2r,\±r^o)± J2e(l±e)cos$;,(2r?, + 4)]. (11) 
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The results for Tqq are consistent with the predictions of refs. [19, 2C, 21], with the 



possible exception of model [^, which at large Q is somewhat higher than the 
data. 

The SCHC-breaking effects observed are not large. This can be seen, for instance, from 
the fact that the SCHC-based relations r\_i = — Im{r^_]^}, Re{rfQ} = — Im{r^Q} are 
satisfied within errors: 

r|_i + Im{r2_J = 0.003 ± 0.016 (stat.) ± 0.030 (syst.) (BPC), 
r|_i + Im{r2_J = -0.037 ± 0.028 (stat.) ± 0.038 (syst.) (DIS), 

Reirl^} + lra{r%} = 0.001 ± 0.005 (stat.) ± 0.013 (syst.) (BPC), 
M^io} + Im{r^o} = 0.006 ± 0.012 (stat.) ± 0.011 (syst.) (DIS). 

The relation 1 — Tqq — 2rl_-^ = (which requires NPE in addition to SCHC) is also 
satisfied, but only within approximately 1.6 standard deviations (summing statistical and 
systematic uncertainties in quadrature) for the DIS data: 



1 - - 2r|„i = 0.060 ± 0.021 (stat.) ± 0.047 (syst.) (BPC), 
l-r^oo- 2r|_i = 0.079 ± 0.032 (stat.) ± 0.036 (syst.) (DIS). 



The ratio of the electroproduction cross section for longitudinal photons to that for 
transverse photons, R, was also determined under a much weaker set of assumptions 
than SCHC (Eqs. (0)-(D). The BPC data give R = 0.368l[i:|]2o (stat.) to^oH (syst.), 
to be compared with the SCHC-based value R = 0.377 ± 0.021 (stat.) ± 0.021 (syst.). 
The DIS data yield R = 2.58l[j;{^ (stat.) toU (syst.) while the SCHC-based value is 
R = 2.661*^0^5 (stat.) ta2i(syst.). The new results differ from the value derived from 
the SCHC hypothesis by less than 3% when averaged over either the BPC or DIS data 
samples. It can therefore be concluded that the results for R presented in Ref . are not 
significantly influenced by the SCHC assumption used in obtaining them. 

The size of the SCHC-breaking effects can also be quantified by evaluating the ratios of 
the helicity-single-fiip and helicity-double-fiip amplitudes to the helicity-conserving am- 
plitudes: 

• The ratio of Tqi (for the production of longitudinally polarised mesons from 
transverse photons) to the helicity-conserving amplitudes can be estimated as (cf. 
Eq. (I)) 

roi = , - (12) 
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which holds under the same assumptions specified for Eq. (^. Equation (\L2^ gives 
Toi = (6.9±1.3 (stat.)±2.4 (syst.))% for the BPC data and tqi = (7.9 ±1.6 (stat.) ± 
2.0 (syst.))% for the DIS data. The HI result for this quantity is (8 ± 3)% |ig. 

Likewise, the ratio 

r,, = l^iol ^ \Re{rro} + Re{rl,}\ 

'■TooP + |Trl2 ^ 



'001 -r I -til I 

can be evaluated, where Tio is the amplitude for the production of transversely 
polarised mesons from longitudinal photons. Equation ( ]T3D is exact if NPE 
holds and Too and Tio have the same phase. The value tiq = (2.9 ±1.4 (stat.) ± 
3.0 (syst.))% is found for the BPC sample and m = (2.7±1.9 (stat.) ±3.1 (syst.))% 
for the DIS data. These results are consistent with the expectation of a severe 
suppression of the production of transversely polarised p° mesons from longitudinal 
photons p^ . 



The size of the helicity-double-fiip amplitude Ti_i can be derived assuming NPE 
and the same phase for Ti_i and Tn: 

n-i = , - (14) 

|Too|2±|rn|2 ./2^ 



which gives ri_i = (4.8 ±1.1 (stat.) ± 2.6 (syst.))% in the BPC analysis and 
ri_i = (1.4 ± 2.7 (stat.) ± 5.9 (syst.))% for the DIS sample. 



Finally, the data can be used to verify the validity of the NPE hypothesis. The left- 
hand-side of Eq. (|Tg) is 0.058 ± 0.027 (stat.) ± 0.074 (syst.) for the BPC sample and 
0.091 ± 0.043 (stat.) ± 0.065 (syst.) for the DIS sample. Both results are consistent with 
the NPE expectation of zero. 



9 Summary and conclusions 



Exclusive electroproduction of mesons has been measured at HERA in two ranges, 
0.25 < Q'^ < 0.85 GeV^ and 3 < < 30 GeV^, and the angular distributions of the decay 
pions have been studied. The low-Q^ data span the range 20 < < 90 GeV; the high-Q^ 
data cover the 40 < W < 120 GeV interval. Both samples extend up to \t\ = 0.6 GeV^. 
The available statistics enable a determination of the 15 combinations of spin-density 
matrix elements, and r^, and allow a check of the extent to which s-channel helicity 
conservation holds at HERA energies for this process. An overall comparison of the data 
with the SCHC hypothesis leads to a very poor x^, /ndf = 81.1/12, for the BPC sample 
and 'X^ /ndf = 79.2/12 for the DIS sample. The deviation from SCHC is most evident in 
the $/j distribution. 
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In the low-Q sample, five of tlie r?^ and r^^ combinations deviate from the null value 
expected if SCHC is assumed. The combination Tqq is measured to be Tqq = 0.051 ± 
0.010 (stat.) ± 0.018 (syst.), which indicates the production of longitudinally polarised 
(i.e. helicity zero) p° mesons from transverse photons. In addition, Rejr^'o} = 0.034 ± 
0.006 (stat.) ±0.009 (syst.), = -0.040 ±0.009 (stat.) ± 0.019 (syst.), = -0.039 ± 
0.009 (stat.) ± 0.021 (syst.) and Imjr^o} = 0.039 ± 0.008 (stat.) ± 0.012 (syst.). The 
measured values for Re{r°o} and r^ti are consistent with the results of earlier fixed-target 
experiments. Some variation of these results with M^^r is observed, possibly indicating 
a contribution from non-resonant tt^tt^ production and its interference with resonant p° 
production. 

In the high-Q^ sample, the only coefficient significantly different from the SCHC expec- 
tation is Tqq, for which the value r^Q = 0.095 ± 0.019 (stat.) ± 0.024 (syst.) is measured. 

00 



This value of Tqq is consistent with the predictions of the pQCD based models |jT9|, ^H, pT 



The results for all coefficients r^^ and at large Q are in good agreement with those 



found by HI ^ 



The ratio R of the cross sections for p° production from longitudinal and transverse 
photons was also determined without assuming SCHC. The results thus found differ from 
those derived from the SCHC hypothesis by less than 3% when averaged over either the 
BPC or DIS data samples. The results for R presented in Ref. [§] are thus not significantly 
influenced by the SCHC assumption used in obtaining them. 

Finally, the data support the hypothesis of natural parity exchange in p° electroproduc- 
tion. 

In conclusion, a small breaking of SCHC, less than 10% in the amplitudes, is a charac- 
teristic feature of exclusive meson production also at the large values of W accessible 
at HERA. In the high-Q^ region, the effect is quantitatively reproduced by perturbative 
QCD calculations. 
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^*^l'10J 2 Nr+eNL ^ Nt+^-Nl ^2 Nt+^Nl 

„04 _ e|Tiop , Re{riiT*_J 
' 1-1 Nr+eNL ' Nr+tNL 



1 _ Re{Ti-iT*J 
' 11 Nr+eNi^ 

IToiP 



^00 Nr+eNL 

1 Re{TnT*_ J _^ i Re{ri-ir*J 



Re-fri I = iilflliliozi!: + 

_ 1 ITiiP I 1 |Ti_ip 
' 1-1 2 Nr+eNL ~^ 2 Nr+eNL 

J^r^2 I _ 1 Re{rnT*_ J i Re{Ti-iT*} 
■^^^^l' 10/ 2 Nr+eNL ~'~ 2 Nr+eNL 

T™r^2 1 _ 1 l^iiP I 1 |Ti-ip 
■^^^^I'l-lJ 2 Nr+eNL ~ 2 Nr+eNL 

5 ^ 1 Re{Tior*} _ 1 Rc{riuT*_J 
'11 ^ Nr+eNL V2 Nr+eNL 

5 ^ V2Re{Toor*J 
' 00 Nr+eNL 

RpJrS 1 = 1 Re{Tiir*„} i Re{Tiur*} l Re{Ti_iT*J 

J^*=\'10/ 2V2 Nr+eNL V2 Nr+eNL 2V2 Nr+eNL 

_ I Rc{TiiT*i(,} 1 Re{TioTliJ 
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1-1 v/2 iVT+eA^z, ^ \/2 Nr+eNL 

1 Re{TiiT*„} 1 Re{Ti-ir*„} 



M„6 \ 1 Kej-iii-inni 1 H-e|Ji-Uni 
'10/ 2V2 Nr+eNL 2V2 Nr+eNL 

T™r„6 -._ 1 Re{T-ior*J 1 Re{rior*„} 

^"n'l-l/ ^ Nr+eNL ^2 Nr+eNL 

where TV^ = iTooP + 2|rioP TVt = iTnP + iT^iP + |T( 
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Table 1: The 15 combinations of spin-density matrix elements r°^, r^, expressed in terms 
of the helicity amplitudes; natural parity exchange in the t channel is assumed for r^ti, 
r^o, rl_^, Ini{r2_J, Refr^J, Im{r6J. 
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ZEi;S 1995 (BPC) 


ZEi;S 1995 (DIS) 


' 00 


0.272 ± 0.011 ± 0.011 


0.725 ± 0.012 ± 0.017 




034 ± 006 ± 009 


0.013 ± 0.010 ± 0.022 


04 
' 1-1 


-0 040 ± 009 ± 019 


000 ± Oil ± 008 


' 11 


-0 039 ± 009 ± 021 


-0.006 ± 0.012 ± 0.026 


^00 


004 ± 015 ± 038 


-0 013 ± 041 ± 076 


-^^^1' 10/ 


-0 019 ± 008 ± 013 


-0 036 ± 015 ± 015 


1—1 


0.334 ± 0.011 ± 0.023 


0.098 ± 0.016 ± 0.016 


Im{r2j 


0.039 ± 0.008 ± 0.012 


0.008 ± 0.014 ± 0.031 




-0.331 ± 0.011 ± 0.020 


-0.135 ± 0.017 ± 0.035 




0.008 ± 0.004 ± 0.012 


0.018 ± 0.005 ± 0.012 


' 00 


0.051 ± 0.010 ± 0.018 


0.095 ± 0.019 ± 0.024 


Re{rfo} 


0.142 ± 0.004 ± 0.012 


0.142 ± 0.007 ± 0.008 


rti 


-0.010 ± 0.006 ± 0.008 


-0.003 ± 0.008 ± 0.005 


Im{r?o} 


-0.141 ± 0.003 ± 0.005 


-0.136 ± 0.007 ± 0.008 


Im{,tJ 


0.01 1 ± 0.006 ± 0.007 


0.009 ± 0.008 ± 0.018 



Table 2: The 15 combinations of spin-density matrix elements r°^, r^^, as obtained from 

the BPC and DIS data sets. Statistical and systematic uncertainties are given separately. 
The BPC data cover the kinematic range 0.25 < < 0.85 GeV^ {{Q^) = 0.41 GeV^), 
20 < < 90 GeV {{W) = 45 GeV), 0.6 < M^^ < 1.0 GeV and \t\ < 0.6 GeV^ 
{{\t\) = 0.14 GeV^). The DIS data cover the kinematic range 3 < Q'^ < 30 GeV^ 
1{Q^) = 6.3 GeV2), AO <W < 120 GeV {{W) = 73 GeV), 0.6 < < 1.0 GeV and 
\t\ < 0.6 GeV2 {{\t\) = 0.17 GeV^). 
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ZEUS 1995 (BPC) 




20 < < 45 GeV 


45 < < 90 GeV 


'00 


0.293 ± 0.014 ± 0.016 


0.257 ± 0.017 ± 0.021 


Rc{r%} 


0.034 ± 0.009 ± 0.008 


0.027 ± 0.010 ± 0.005 


' 1-1 


-0.032 ± 0.011 ± 0.009 


-0.056 ± 0.013 ± 0.020 


1 


-0.040 ± 0.012 ± 0.006 


-0.038 ± 0.015 ± 0.007 


^00 


0.004 ± 0.021 ± 0.030 


0.003 ± 0.025 ± 0.051 


r 1 1 

Re{riJ 


—0.024 ± 0.011 ± 0.016 


-0.012 ± 0.013 ± 0.005 


ri_i 


0.328 ± 0.015 ± 0.014 


0.345 ± 0.018 ± 0.018 


Im{r2J 


0.037 ± 0.010 ± 0.013 


0.038 ± 0.012 ± 0.006 


Im{r?_i} 


-0.313 ± 0.015 ± 0.023 


-0.348 ± 0.018 ± 0.022 




0.015 ± 0.006 ± 0.008 


0.001 ± 0.007 ± 0.010 


'00 


0.039 ± 0.013 ± 0.016 


0.050 ± 0.016 ± 0.007 


Re{rfo} 


0.143 ± 0.005 ± 0.008 


0.142 ± 0.006 ± 0.013 




-0.013 ± 0.008 ± 0.007 


-0.008 ± 0.009 ± 0.010 


Im{^?o} 


-0.150 ± 0.004 ± 0.004 


-0.131 ± 0.006 ± 0.007 


lm{ri_i} 


0.011 ± 0.007 ± 0.004 


0.022 ± 0.009 ± 0.007 



Table 3: The 15 combinations of spin-density matrix elements r°^, r^"^, as obtained from 
the BPC data in two W intervals. Statistical and systematic uncertainties are given 
separately. The data cover the kinematic range 0.25 < < 0.85 GeV^, 20 < 1^ < 
90 GeV, 0.6 < M^^ < 1.0 GeV and \t\ < 0.6 GcV^. The average W values for the two 
intervals are {W) — 31 GeV and {W) = 61 GeV, respectively. 
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ZEUS 1995 (BPC) 




0.25 < Q < 0.4 GeV^ 


0.4 < Q < 0.85 GeV^ 


'00 


0.228 ± 0.016 ± 0.022 


0.324 ± 0.015 ± 0.012 


Rc{r%} 


0.025 ± 0.009 ± 0.004 


0.040 ± 0.010 ± 0.009 


' 1-1 


-0.053 ± 0.012 ± 0.003 


-0.028 ± 0.012 ± 0.006 


1 


-0.039 ± 0.014 ± 0.018 


-0.033 ± 0.013 ± 0.009 


^00 


-0.031 ± 0.022 ± 0.022 


0.038 ± 0.023 ± 0.045 


r 1 1 

Re{riJ 


-0.020 ± 0.011 ± 0.007 


-0.016 ± 0.012 ± 0.015 


ri_i 


0.347 ± 0.016 ± 0.013 


0.329 ± 0.016 ± 0.012 


Im{r2J 


0.035 ± 0.010 ± 0.005 


0.045 ± 0.011 ± 0.008 


Im{r?_i} 


-0.353 ± 0.016 ± 0.018 


-0.302 ± 0.016 ± 0.031 




0.009 ± 0.006 ± 0.010 


0.007 ± 0.006 ± 0.004 


'00 


0.042 ± 0.014 ± 0.011 


0.055 ± 0.015 ± 0.009 


Re{rfo} 


0.139 ± 0.005 ± 0.012 


0.149 ± 0.006 ± 0.008 




-0.006 ± 0.008 ± 0.007 


-0.011 ± 0.008 ± 0.009 


Im{^?o} 


-0.128 ± 0.005 ± 0.006 


-0.159 ± 0.005 ± 0.005 


lm{ri_i} 


0.007 ± 0.008 ± 0.006 


0.025 ± 0.008 ± 0.008 



Table 4: The 15 combinations of spin-density matrix elements r°^, r^^, as obtained from 
the BPC data in two Q'^ intervals. Statistical and systematic uncertainties are given 
separately. The data cover the kinematic range 0.25 < < 0.85 GeV^, 20 < 1^ < 
90 GeV, 0.6 < M^^ < 1.0 GeV and \t\ < 0.6 GeV^. The average values for the two 
intervals are (Q^) = 0.32 GeV^ and (Q^) = 0.54 GeV^, respectively. 
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ZEUS 1995 (BPC) 




U < c < U.i Kje\ 


U.i < \t\ < U.zo KjieW 


U.zo < t < U.D IjeV 


' 00 


no<^i -L nni/i -L nni7 
U.zDi ± U.Ui4 ± U.Ui( 


n 07Q -L n non -l n ni ft 
U.Zio ± U.UzU ± U.UiD 


n QQO -L n non -l n no7 
u.ooy ± u.uzy ± u.Uz( 




u.uzy ± u.uuy ± U.UiZ 


n n/io -I- n mo -1- n nnfi 
U.U4Z it U.UiZ it U.UUo 


n n/17 -I- n ni 7 -1- n no/i 
U.U4/ it (j.Vii ± U.Uz4 


^04 
^1-1 


— U.Uoo ± U.Uiz ± U.Ui4 


n nQ/1 1 nm^. i nm" 
— U.Uo4 it U.Uii3 it U.Uio 


nn/is -1- nnin -i- nmo 
— U.U4o it U.Uiy it U.Uiz 


^11 


— U.UZO ± U.Uio ± U.UiO 


— U.U4o ± U.Ui/ ± U.Uiy 


n nt^/i -1- n noo _i_ n noi 
— U.Uy4 it U.Uzz it U.Uzi 


^00 


U.Uzo ± U.UzU ± U.Uo4 


n non _i_ n noo _i_ n npci 
— U.Uzy it U.Uzo it U.Ubi 


n nvi _i_ n n/ip; _i_ n nc^ 
— U.U/i ± U.U4t) ± U.Uoo 




nnio-unnii i n m -i 
— U.Uiz it U.Uii it U.Ui4 


nnQO -4- nni/i -u nnis 
— U.Uoz it U.Ui4 it U.Uio 


n nifi -4- n noi -u n nQO 
— U.Uio it U.UZi it U.Uoz 


^i-i 


0.328 ± 0.016 ± 0.015 


0.333 ± 0.020 ± 0.023 


0.300 ± 0.026 ± 0.020 


Im{^?o} 


0.027 ± 0.010 ± 0.013 


0.038 ± 0.014 ± 0.019 


0.044 ± 0.020 ± 0.030 


Im{r2_,} 


-0.329 ± 0.016 ± 0.025 


-0.319 ± 0.020 ± 0.027 


-0.285 ± 0.030 ± 0.030 


^11 


0.001 ± 0.006 ± 0.008 


0.012 ± 0.008 ± 0.004 


0.021 ± 0.010 ± 0.014 


/i-.5 
'GO 


0.026 ± 0.012 ± 0.025 


0.070 ± 0.020 ± 0.009 


0.117 ± 0.027 ± 0.050 


Re{r5J 


0.137 ± 0.005 ± 0.011 


0.152 ± 0.007 ± 0.012 


0.144 ± 0.010 ± 0.014 


'1-1 


-0.002 ± 0.008 ± 0.008 


-0.012 ± 0.010 ± 0.008 


-0.018 ± 0.014 ± 0.010 


Im{r?o} 


-0.144 ± 0.005 ± 0.009 


-0.135 ± 0.007 ± 0.004 


-0.121 ± 0.011 ± 0.022 


lm{rf_i} 


0.012 ± 0.008 ± 0.008 


0.007 ± 0.010 ± 0.007 


0.009 ± 0.015 ± 0.014 



Tabic 5: The 15 combinations of spin-density matrix elements r^^, r^, as obtained 
from the BPC data in three t intervals. Statistical and systematic uncertainties are given 
separately. The data cover the kinematic range 0.25 < < 0.85 GeV^, 20 < W < 
90 GeV, 0.6 < M^^ < 1.0 GeV and \t\ < 0.6 GeV^. The average |t| values for the three 
intervals are {\t\) = 0.04 GeV^, {\t\) = 0.16 GeV^ and {\t\) = 0.37 GeV^, respectively. 
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ZEUS 1995 (BPC) 




'oo 


Re{rM} 


,„04 
'l-l 






Re{r;} 






Im{rj_j} 


A 


'oo 


Re{4} 






Im{rJ_j} 


'00 


1,000 
































0,079 


1,000 




























.,04 
'l-l 


0,046 


0,054 


1,000 




























0,101 


0,011 


0,584 


1,000 
























J 

'oo 


-0,159 


-0,025 


-0,062 


-0,414 


1,000 






















Mrlo} 


-0,068 


-0,508 


-0,127 


-0,039 


0,043 


1,000 




















rl-i 


-0,327 


-0,201 


-0,082 


-0,082 


0,052 


0,116 


1.000 


















Wn} 


0,105 


0,445 


0,075 


0,062 


-0,039 


-0,293 


-0,095 


1,000 
















Im{rj_i} 


0,312 


0,097 


0,019 


0,009 


-0,029 


-0,069 


0,073 


0,072 


1,000 














r\, 


-0,095 


0,231 


0,014 


-0,075 


0,063 


-0,038 


0,014 


0,011 


0,050 


1,000 












^5 

'oo 


0,248 


0,578 


0,053 


0,087 


-0,223 


-0,452 


-0,150 


0,544 


0,133 


-0,300 


1,000 










Re{rjo} 


0,299 


0,080 


-0,225 


-0,138 


-0,184 


-0,201 


-0,098 


-0,084 


0,023 


-0,004 


-0,044 


1,000 










0,003 


-0,325 


0,052 


0,058 


0,007 


0,117 


-0,015 


-0,012 


-0,068 


-0,479 


-0,051 


0,078 


1,000 






HA} 


-0,308 


0,051 


-0,359 


-0,221 


-0,356 


0,031 


0,009 


-0,056 


-0,057 


0,044 


0,034 


0,056 


-0,054 


1,000 




Im{ri_i} 


-0,009 


0,337 


-0,046 


-0,074 


0,014 


0,062 


-0,072 


0,169 


0,035 


0,443 


0,041 


-0,040 


-0,251 


0,055 


1,000 



Table 6: Correlation matrix for the 15 combinations of spin-density matrix elements r' 
r"i^, as obtained from the BPC data. 
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ZEUS 1995 (DIS) 




'oo 


Re{rM} 


,„04 
'l-l 




'''oo 


Re{r;} 






Im{rj_j} 


A 


'oo 


Re{4} 






Im{rJ_j} 


'oo 


1,000 
































-0,014 


1,000 




























.,04 
'l-l 


0,020 


-0,081 


1,000 




























0,061 


-0,097 


0,273 


1,000 
























J 

'oo 


-0,011 


-0,042 


0,070 


-0,381 


1,000 






















Mrlo} 


0,086 


-0,219 


0,146 


-0,056 


-0,086 


1,000 




















rl-i 


-0,203 


0,059 


-0,060 


-0,050 


0,023 


0,196 


1.000 


















Wn} 


0,017 


0,176 


0,013 


-0,032 


-0,061 


0,152 


0,099 


1,000 
















Im{rj_i} 


0,182 


-0,008 


-0,028 


0,041 


0,005 


0,170 


0,433 


0,095 


1,000 














r\, 


-0,257 


0,292 


-0,050 


-0,174 


0,088 


-0,038 


-0,054 


0,166 


-0,086 


1,000 












^5 

'oo 


0,172 


0,422 


-0,027 


0,067 


-0,221 


-0,333 


-0,004 


0,254 


0,042 


-0,279 


1,000 










Re{rjo} 


-0,131 


0,331 


-0,460 


-0,170 


-0,322 


-0,221 


0,023 


-0,029 


-0,162 


0,054 


-0,004 


1,000 










0,015 


-0,270 


0,309 


-0,022 


-0,013 


0,179 


-0,107 


0,195 


-0,094 


-0,156 


-0,028 


-0,038 


1,000 






HA} 


0,124 


0,076 


-0,432 


-0,220 


-0,339 


-0,134 


-0,143 


-0,214 


0,102 


0,021 


0,092 


0,506 


-0,014 


1,000 




Im{ri_i} 


-0,071 


0,244 


0,109 


-0,043 


0,007 


0,205 


-0,045 


0,284 


-0,245 


0,221 


-0,002 


0,001 


0,286 


0,013 


1,000 



Table 7: Correlation matrix for the 15 combinations of spin-density matrix elements r' 
r^^, as obtained from the DIS data. 
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ZEUS 1995 (BPC) 




0.6 < M^^ < 0.77 GeV 


0.77 < M^^ < 1.0 GeV 


'00 


0.311 ± 0.014 ± 0.013 


0.227 ± 0.018 ± 0.016 


Re{r?on 


0.026 ± 0.009 ± 0.010 


0.032 ± 0.010 ± 0.005 


'1-1 


-0.031 ± 0.011 ± 0.006 


-0.062 ± 0.014 ± 0.010 


'•!i 


-0.030 ± 0.012 ± 0.009 


-0.051 ± 0.015 ± 0.007 


'oo 


0.020 ± 0.019 ± 0.042 


-0.011 ± 0.026 ± 0.041 


Re{rlo} 


-0.007 ± 0.010 ± 0.010 


-0.029 ± 0.013 ± 0.010 




0.321 ± 0.014 ± 0.016 


0.369 ± 0.018 ± 0.027 




0.023 ± 0.010 ± 0.010 


0.064 ± 0.012 ± 0.017 


lm{rf_i} 


-0.302 ± 0.014 ± 0.029 


-0.362 ± 0.019 ± 0.016 




0.015 ± 0.006 ± 0.009 


-0.002 ± 0.007 ± 0.006 


' 00 


0.026 ± 0.013 ± 0.019 


0.077 ± 0.015 ± 0.017 


Re{r^J 


0.152 ± 0.005 ± 0.008 


0.129 ± 0.007 ± 0.012 


rti 


-0.010 ± 0.007 ± 0.007 


-0.009 ± 0.009 ± 0.013 


Im{r^o} 


-0.155 ± 0.004 ± 0.004 


-0.125 ± 0.006 ± 0.006 




0.019 ± 0.007 ± 0.006 


0.007 ± 0.009 ± 0.008 



Tabic 8: The 15 combinations of spin-density matrix elements r°^, r^"., as obtained from 
the BPC data in two M^^^^ intervals. Statistical and systematic uncertainties are given 
separately. The data cover the kinematic range 0.25 < < 0.85 GeV^, 20 < W < 
90 GeV, 0.6 < M^^ < 1.0 GeV and \t\ < 0.6 GeV^. 
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Figure 1: Schematic diagram of the reaction ep — > ep°p, indicating the kinematic variables 
used in this analysis. 
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Figure 2: Schematic diagrams of (a) the process ep — > ep^p in the 7*p centre-of-mass 
system, and (b) the decay of the p° in its rest frame. Three angles suffice to describe the 
reaction: the azimuthal angle between the scattering plane and the production plane, <^y^\ 
and the two p° decay angles, 0/^, the azimuthal angle between the production and decay 
planes, defined in either the 7*]? system or in the rest frame; and 9^, which is the polar 
angle of the positively-charged decay product defined with respect to the direction of the 
p° momentum vector in the 7*p system, or, equivalently, the direction opposite to the 
momentum- vector of the final-state proton in the rest frame of the p° meson. This choice 
of the spin-quantisation axis defines the helicity frame. 
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Figure 3: Acceptance as a function of cosOh, (f)^ and for the BPC and the DIS data 
samples. The bars indicate the statistical uncertainties. 
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Figure 4: Observed distributions for cos 9h (a), 0/i (b), (c) and iph = (ph — (d, 
e, f) of the reconstructed BPC data (points) and the reconstructed Monte Carlo events 
(histograms). The sohd histograms correspond to the Monte Carlo data reweighted with 
Eq. d^), in which the results of the present analysis were used for the 15 combinations 
of matrix elements r^^, r^^. The dashed histograms correspond to the SCHC hypothesis. 
The distributions are not corrected for acceptance. The error bars indicate the statistical 
uncertainty. The statistical uncertainty of the simulated distributions is negligible. 
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Figure 5: Observed distributions for cos6'/i (a), (ph (b), (c) and i'h = (ph — (d, 
e, f) of the reconstructed DIS data (points) and the reconstructed Monte Carlo events 
(histograms). The sohd histograms correspond to the Monte Carlo data reweighted with 
Eq. d^), in which the results of the present analysis were used for the 15 combinations 
of matrix elements r^^, r^^. The dashed histograms correspond to the SCHC hypothesis. 
The distributions are not corrected for acceptance. The error bars indicate the statistical 
uncertainty. The statistical uncertainty of the simulated distributions is negligible. 
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Figure 6: The 15 combinations of spin-density matrix elements r°^, r^^, as obtained from 
the BPC and DIS data sets as a function of . The full symbols indicate the present 
results; the open points indicate the HI results |jTO|. The inner error bars indicate the 
statistical uncertainties, the outer the statistical and systematic uncertainties summed 
in quadrature. The continuous curves are the prediction of the model calculation of 
Ref. the dotted curves the prediction of Ref. [^. The dashed lines indicate the 
SCHC expectation. The predictions of Ref. ||2T[ and of SCHC coincide for r^l^ and r^. 
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Figure 7: The 15 combinations of spin-density matrix elements r°^, r^, as obtained from 
the BPC data as a function of t. The inner error bars indicate the statistical uncertainties, 
the outer bars the statistical and systematic uncertainties summed in quadrature. The 
horizontal bars indicate the size of the bins. The dashed lines indicate the prediction of 
SCHC, where available. The data cover the kinematic range 0.25 < < 0.85 GeV^, 
20 < < 90 GeV, 0.6 < M^^ < 1.0 GeV and \t\ < 0.6 GeV^. 
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Figure 8: The 15 combinations of spin-density matrix elements r?^, r^^, as obtained 



from the DIS data compared with the predictions of the model calculation of Ref. |119 
(continuous histogram) and of SCHC (dashed histogram, see text for details). The solid 
points indicate the present results; the open points indicate the HI data [|10|. The inner 
error bars indicate the statistical uncertainties, the outer the statistical and systematic 
uncertainties summed in quadrature. The ZEUS data cover the kinematic range 3 < < 
30 GeV^ ((g2) = 6.3 GeV^), AO < W < 120 GeV {{W) = 73 GeV), 0.6 < M^^ < 1.0 GeV 
and \t\ < 0.6 GeV^. The HI data cover the range 1 < Q'^ < 60 GeV^, 30 < W < 140 GeV 
and \t\ < 0.5 GeV^. The prediction of Ref. [0 is for = 10 GeV^ and W = 100 GeV. 
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